Non-linear resonance in the accretion disk of a millisecond pulsar 
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ABSTRACT 

Twin quasi-periodic millisecond modulations of the X-ray flux (kHz QPOs) have 
recently been reported from an accreting 2.5 ms X-ray pulsar. We identify modes 
of disk oscillations whose frequencies are in agreement with the observed ones when 
the rotating neutron star is modeled with realistic equations of state. The frequency 
difference of the twin QPOs, equal to about one half of the neutron-star spin rate, 
clearly indicates that resonant oscillations of the accretion disk have been observed. 
Similar non-linear resonances may also be spontaneously excited in the accretion disk. 
The two QPO frequencies in the pulsar system are close to a 5:7 ratio and this suggests 
a link with the QPOs in black hole systems, where frequency ratios of 2:3 and 3:5 have 
been reported. 

Subject headings: Stars: neutron - X-rays: general 



1. Introduction 



A millisecond pulsar that accretes matter from a binary stellar companion provides a unique 
probe of the physics of accretion flow a few Schwarzschild radii away from the neutron star. In 
contrast with the case of the slowly rotating steady X-ray pulsars, the magnetic field of the (X-ray 
transient) millisecond pulsar is too weak to prevent the accreting gas from orbiting the neutron 
star very close to its surface. Hence, in such a pulsar oscillations in the inner parts of the disk may 
be observed at their characteristic frequency of several hundred Hertz. The pulsar is expected to 
disturb the accretion disk at its spin frequency. That it does, is demonstrated by the discovery in 
the 2.5 ms accreting pulsar of a frequency difference between the two QPOs equal to one-half the 
pulsar spin rate. 
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We point out that the two QPO frequencies themselves may be commensurable with the spin 
frequency and with each other. The two simultaneously observed QPO frequencies are close to 
500 Hz and 700 Hz, i.e., they may be in a 5:7 ratio. Twin QPOs have been proposed to be a 
manifestation of non-linear resonance in relativistic accretion disks (Kluzniak and Abramowicz 
2001), and the resulting suggestion that twin QPO frequencies should be in the ratio of small 
integers has subsequently been confirmed for a number of black hole transients (Abramowicz and 
Kluzniak 2001; Remillard et al. 2002; Kluzniak and Abramowicz 2002; McClintock and Remillard 
2003) and for the bright steady source Sco X-1 (Abramowicz et al. 2003a). 

The two QPOs discovered in the 2.5 ms accreting pulsar (Wijnands et al. 2003) are similar to 
QPOs previously detected in many unpulsed X-ray sources, in which an accretion disk is thought to 
feed matter and angular momentum to the neutron star or black hole (van der Klis et al. 2000; Lewin 
et al. 1995). Such QPOs in black hole systems were intrerpreted as oscillations of the accretion disk 
(Wagoner 1999; Kato 2001), but a direct proof of this was lacking, and for neutron star systems 
other models had been proposed (van der Klis et al. 2000). SAX J1808.4-3658 is the first kHz QPO 
source in which the rotational period of the neutron star has been unambiguously determined, and 
the first in which the frequency difference of the twin QPOs is known to be a subharmonic of the 
stellar rotation rate. For the first time, we have direct evidence that kHz QPOs are caused by disk 
oscillations. 



2. Non- linear resonance and epicyclic frequencies 

The presence of a sub-harmonic frequency is a clear signature of non-linear resonance (Nayfeh 
and Mook 1979). Theodor von Karman observed in 1940 that high-frequency vibrations of an 
airplane engine excite lower-frequency resonances in the airframe (von Karman 1940; Nayfeh and 
Mook 1979). A tragic example is known when a subharmonic resonance was excited in an airplane 
wing, which in turn excited a subharmonic resonance in the rudder, at 1/2 the wing eigenfrequency 
(Lcfschctz 1956; Nayfeh and Mook 1979). In the present case, the 2.5 ms pulsar SAX J1808. 4-3658 
plays the role of the engine, while the accretion disk (whose response is observed as the quasi- 
periodic oscillations) is the airframe with its own set of known eigenfrequencies. The fact that the 
QPO frequency difference coincides with one half the known pulsar spin frequency, 

(694 ± 4) Hz - (499 ± 4) Hz = (195 ± 6) Hz « ^(401 Hz), 

follows directly from the nature of the system. The structure of the accretion disk is determined 
by the non-linear equations of hydrodynamics, and the effective gravitational potential (in which 

the accreting gas moves) is not harmonic, so a non-linear response is expected. Below, we discuss 
in detail how this comes about for an accretion disk in the space-time of a rotating neutron-star. 

The motion of a test particle in nearly circular orbits close to the equatorial plane can be 
decomposed into three components, circular planar motion at the orbital frequency SI = 27ri/orb) 



-3- 



harmonic radial motion at the radial epicyclic frequency n = and harmonic vertical motion 

at the meridional epicyclic frequency C, = 27rzvvcrt- In Newtonian theory of spherically symmetric 
gravitating bodies the three frequencies coincide, but in Einstein's gravity k <n. (The relativistic 
precession of the perihelion of Mercury occurs at the rate Vl — n.) The same three frequencies are 
important in a discussion of fluid motion about a gravitating body in general relativity (Kato 1998). 
The accretion disk is a body of hot gas that is supported against infall primarily by rotation and 
is nearly in hydrostatic equilibrium (Lewin et al. 1995; Kato 1998). Theories of accretion (Shakura 
and Sunyaev 1973; Jaroszyhski et al. 1980) admit solutions in which the disk thickness is much 
smaller than its radial extent, but also solutions in which the disk has the geometry of a torus. 
Like other extended bodies in equilibrium, the disk is capable of motion in a variety of modes. In 
the linear regime, small oscillations of geometrically thin accretion disks, as well as waves in this 
body of fluid, have been extensively studied in the Kerr metric (Wagoner 1999; Kato 2001). The 
radial vibrations of two-dimensional models of geometrically thick disks (accretion torii) have been 
investigated numerically in the Schwarzschild metric (Rczzolla et al. 2003). Qualitatively, these 
black-hole metrics are similar to the metrics of neutron stars. 

Wc have studied numerically the response to a transient external perturbation of an ideal-gas 
torus that is initially in equilibrium rotation about a (Schwarzschild) black hole or neutron star. 
We find that the torus performs harmonic oscillations both in the radial and vertical directions, 
even if the impulse imparted to it at the beginning of the computation is purely radial. Further, 
the vertical oscillations have variable amplitude. Both these effects speak of a non-linearity in 
the system. The vertical and radial displacements of the center of the torus (defined as the point 
of maximum pressure in a meridional cross-section) are shown as a function of time in Fig. 1. 
The radial displacement has been rescaled to fit the graph — in the simulation described here, the 
amplitude of radial motion is actually larger than the amplitude of vertical motions. 

For the numerical simulation (Fig. 1), we use a Newtonian 2-d smooth particle hydrodynamics 
(Monaghan 1992, SPH) code (Lee and Ramirez-Ruiz 2002) simulating ideal gas with adiabatic index 
4/3 in a pseudo-potential (Kluzniak and Lee 2002) ^KL{r) = [1 — exp(rms/'')]G'M/rms, which 
reproduces the Schwarzschild ratios of the orbital and epicyclic frequencies, (/Q, = 1, (k/C)^ = 
l—rms/f- Here, r„is = 6M(G/c^) is the radius of the marginally stable orbit. A torus in equilibrium, 
with ro = 12.25M(G/c^), was perturbed at time t = by imparting to it a radial velocity field 
with magnitude proportional to y^rms/r. The frequencies of the resulting motions are found to 
be inversely proportional to the central mass M, whose value was adjusted to match the observed 
QPO frequencies. 

The two frequencies observed in this numerical simulation coincide, within errors, with the 
known epicyclic frequencies, k and (, at the equilibrium radial position, ro, of the locus of maximum 
pressure, rc- We conclude that a slender torus responds to a radial external perturbation with 
oscillations occuring at frequencies equal to the two epicyclic frequencies, ^(ro) and CC'^o); at that 
radius where the pressure of the torus in equilibrium is highest. An important implication is that 
these frequencies can vary if the torus varies with time as it accumulates less or more of the matter 
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fiowing through it. 

For the radial perturbation appHed in the simulation, the amplitude of vertical motion is 
smaller than that of the radial motion. Presumably the ratio of the amplitudes depends on the 
perturbation applied. Quantifying this requires further study. A tilted magnetic dipole may perturb 
the disk quite strongly in the vertical direction. In the remainder of this section we assume that the 
perturbation is such that only one of the two motions is typically excited with a large amplitude, 
and that the other would not be as easily discernible in an X-ray observation, unless its initially 
much smaller amplitude were amplified by a resonance. 

A pulsar is a rotating neutron star with a strong magnetic dipole not aligned with the rotation 
axis. As a result, in the 2.5 ms accreting pulsar (Wijnands and van der Klis 1998) the accretion 
disk suffers a periodic disturbance at the spin frequency ui = 401 Hz. Observations indicate the 
presence of a QPO varying in frequency between 280 and 750 Hz (Wijnands et al. 2003). We 
interpret this as one of the two epicyclic frequencies, uq, and expect that a second QPO will be 
present at the second epicyclic frequency, when a resonance occurs at frequency v^es between the 
corresponding disk oscillation and the disturbing forces. 

Twin QPOs have been observed in this source on only one occasion, and it is not clear to us 
whether the single variable frequency i^o present in the remaining observations corresponds to the 
upper or the lower QPO, i.e., to Uyert or to Ur (but see below). A simple and general reason why 
the difference is equal to one half the spin frequency, v^es — J^o = J^i/2, and not the spin frequency, 
vi, can be given if i/o = u,-- Then i^j-cs = invert at & particular radius, and the resonance would occur 
in an oscillator without quadratic anharmonicity terms. The restoring force in the vertical motion 
of the torus is obtained by taking the z derivative of the effective potential expanded about the 
equatorial plane. Because the potential is symmetric under reflection in this plane, this expansion 
has only even powers of the cylindrical co-ordinate z, V{z) = Vq + {Cz)'^/2 + /3z'^/4+... The vertical 
motion of the torus can then be described as that of an oscillator with (angular) eigen-frequency 
^, and a cubic anharmonicity in leading order, (3z^. 

Consider resonance in an accretion disk which can be described as a damped oscillator with a 
cubic anharmonicity. In view of the very strong disturbance by the pulsar, the resonance may be 
driven by a combination of the pulsar spin frequency vi and another, much weaker, disturbance 
at frequency 1^2, which we take to be a harmonic (i.e., an integer multiple) of the QPO frequency 
fQ. In the present case the frequencies of the two harmonic driving forces satisfy 1^2 > z^i- The 
oscillator may be in resonance with both of these forces at the same time (Nayfeh and Mook 1979, 
simultaneous resonance) in one of three cases: 

a) Pres ~ 2r^i ± 1/2, or 

b) f res ~ 2z/2 ±1^1, OTC 
C) l^ves ~ (z^2 ± Z^l)/2. 

The single occasion when two QPOs were detected in the 2.5 ms pulsar corresponds to case c), with 
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1^2 = 2i/o, i.e., Vres ~ 1^0 ^ 1^1 Since for neutron stars in general relativity invert > i^r, the upper 

sign corresponds to the case z/j-cs = J^vcit, ^'o = ^^r, the lower sign to the case z/jcs = i^r, ^'o = J^vcit- 
Indeed, the two observed QPO frequencies, are in this relation to the spin frequency of 401 Hz, 
fupper - t'lowcr = 694 Hz - 499 Hz « (401Hz)/2 = iyi/2. 

We note that the two QPOs cannot be separated by the spin frequency in this scheme. A 
separation between the two QPO frequencies of twice the spin frequency, {v^es — = 802 Hz, case 
a with uq = 1/2) is in principle also possible in this scheme if resonance occurs at the very inner 
edge of the disk, at frequencies of about 650 Hz and 1450 Hz, but it does not seem likely that the 
bulk of the torus could reside so close to the neutron star. 

We have numerically computed realistic models of rotating neutron stars and their exterior 
metrics. Fig. 2 presents a typical result consistent with observations of the two QPOs in the 2.5 
ms pulsar. The vertical and radial epicyclic frequencies for a 1.22 solar mass neutron star, modeled 
with equation of state FPS, and rotating at 401 Hz arc illustrated as a function of the radius. The 
case c) resonance discussed above occurs when pressure in the torus has a maximum at radius 
tq = 11.27M(G/c^) = 20.6km. The stellar mass and resonance radius are only indicative of the 
true values for this pulsar, as the frequencies for a flattened torus of large radial extent may differ 
from /«(ro) and C(^o) by several per cent (Rezzolla et al. 2003). 

Wijnands et al. (2003) give phenomenological arguments in favor of the frequency being 
the upper QPO. In this case, in addition to the resonance described below, a resonance at the 
combination frequency Vj- = i/'yert ~ could also occur in principle if the radial oscillator with its 
quadratic anharmonicity were perturturbed simultaneously at the spin frequency and the frequency 
of the vertical oscillation. This would happen at ro = 16 km for the same model of the neutron star, 
when the upper QPO frequency has the (as yet unobserved value) 1015 Hz and the lower frequency 
is 615 Hz (dashed vertical line in Fig. 2). 

To compute the frequencies in Figs. 2 and 3, we constructed numerical models of rotating 
neutron stars and their exterior metrics using the relativistic code of Stergioulas and Friedman 
( 1995), which solves Einstein's field equations for arbitrarily large rotation rate in an integral 
form. Details of the numerical method and extensive accuracy tests can be found in (Stergioulas 
2003). The results presented here are not very sensitive to the choice of equation of state (EOS) of 
neutron-star matter. 

3. Parametric resonance 

So far we have been discussing a resonance directly caused by an external disturbance. But 
when the frequencies n{r) and C(r) happen to be in particular ratios, large-amplitude motions of 
the torus may be excited spontaneously. This mechanism, and specifically parametric resonance 
between the two modes of oscillation, has been suggested (Kluzniak and Abramowicz 2002) as the 
origin of the high-frequency QPOs in black hole systems, where the frequency ratios 5:3 and 3:2 
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have been noted (Abramowicz and Kluzniak 2001; Remillard et al. 2002; Kluzniak and Abramowicz 
2002; McClintock and Remillard 2003), as well as in the candidate neutron-star system Sco X-1 
(Abramowicz et al. 2003a, b). The presence of subharmonics in some of the black hole systems has 
also been discussed (Remillard et al. 2002; Kluzniak and Abramowicz 2003). On the other hand, 
in several X-ray bursters the frequency difference between the twin QPOs has been reported to be 
close to the suspected spin frequency (or one half of it) of the neutron star (van der Klis et al. 2000). 
It is interesting to note that in addition to being separated by one-half the pulsar spin frequency, 
the two QPO frequencies in SAX J1808. 4-3658 are approximately in a 7:5=1.4 ratio. The system 
may thus be the missing link between the QPOs in black hole systems and X-ray bursters. It is 
known that if the frequency difference of two modes is accidentally close to that of an external 
disturbance, a resonance involving the three frequencies may occur (Kato 1974). In neutron star 
systems this could favor the appearance of strong QPOs with a frequency difference selected by the 
perturbing neutron-star rotation, while in black hole disks only the spontaneous resonance would 
occur. 

Parametric resonance between two oscillations is possible in accretion disks because the cou- 
pling of modes is non-linear in hydrodynamics. For thin disks, this can lead to excitation of two 
modes whose frequency ratio is closer to 1.4 than to 1.5. In the Schwarschild metric, the damping 
or excitation of ^f-modes with azimuthal dependence ex.p(im(p) has been considered by Kato (Kato 
2003), who found that an m = mode and an m = 1 mode will undergo resonant amplification 
when interacting with a non-rotating one-armed (m = 1) stationary warp, assumed to be present 
in the accretion disk. We assume that the same modes are unstable also in a gravitational field 
that is not spherically symmetric, and compute their frequencies for the metric of rotating neutron 
stars. The sum of the two frequencies turns out to be i^orh, and if they arc identified with the two 
observed QPO frequencies this constrains the mass of the 2.5 ms pulsar to be less than 1.5 solar 
mass. The two mode frequencies, z^ower and i^upper, as well as z/orb) ^r, and i^orb ~ are shown in 
Fig. 3 for a 1.49 solar mass neutron star, spinning at a frequency of 406 Hz. The resonance occurs 
in the region where J/iowcr < '^r , with the largest modulation of X-ray flux expected when the two 
frequencies are evaluated at the edge of the resonance zone, i.e., when they are z^ower ~ and 

Z^upper ~ Z^orb ^r- 

The two unstable modes in Fig. 3 were taken, following (Kato 2003), to be m = 1, n = 1, and 
m = 0, n = 1, g-modes excited by parametric resonance, which is mediated by a warp assumed 
to be present in the thin accretion disk. For a non-rotating warp their frequencies (Kato 2003) 
would be 2forb ~ V^^^vert and \/2i/yert — forb) respectively, with their ratio tending to \/2 = 1.41 
in the limit of non-rotating stars, or in the limit of large radial distance from the star, r — ^ oo. 
The warp can be identified with a c-mode of disk oscillations and is expected to rotate at the 
Lense-Thirring angular frequency 27r(forb — '^vcrt), evaluated at a certain characteristic radius, 
(Silbcrglcit et al. 2001). The two curves plotted in Fig. 3 arc /^iowcr(^) = (v^ — l)i^vcrt('') and 
i'uppev{r) = forb + (1 ~ V^)i'vertif) ■ The actual mode frequencies coincide with these at r = r^, 
and at arbitrary radius differ from these by [I'oThi'r'w) — i^vertifw)] — Woihif) — ^vert{f)], i-e-, by a few 
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percent. For the present purposes, we take the quantity r^, to be an unknown parameter. 

4. Conclusions 

In conclusion, we note again that it is not yet known whether the accretion disk resembles a 
torus or, to the contrary, is geometrically thin. The 2.5 ms pulsar is the first source in which the 
frequencies of the twin QPOs are at once clearly related to the neutron-star spin frequency and are 
in a 1.4 ratio. With one observation of twin frequencies in only one source with a known rotation 
rate (but an unknown mass), it is impossible to determine which of the two pairs of oscillation 
modes discussed here is responsible for the non-linear resonance. The ambiguity will be removed if 
twin QPOs are observed in one or more of the four other accreting pulsars (Markwardt and Swank 
2003) known, where the spin periods are different. 

It is worth noting that the fact that the frequency difference between the two observed kHz 
QPOs is equal to one half the spin frequency of the neutron star is the first clear indication that 
the commonly observed millisecond modulations of the X-ray flux in low mass X-ray binaries are 
caused by oscillatory motions of the accretion disk around the neutron star (or black hole in other 
systems). 
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Fig. 1. — The vertical and radial motions of an accretion torus in response to an impulsive radial 
perturbation at time t = 0. The vertical displacement of the torus (dashed line) about the equatorial 
plane and the radial displacement (continuous line) about the equilibrium radius, rg, of the circle 
of maximum pressure are plotted as a function of time. In this pseudo-Schwarzschild simulation, 
ro = 12.25M(G/c^). It is clear that the torus performs two harmonic motions. We find that their 
frequencies coincide with the meridional and radial epicyclic frequencies, Vyert smd Ur, here 700 Hz 
and 500 Hz, respectively. The ratio of the frequencies is a function of ro, their value a function of 
the mass M of the neutron star. See the text for details. 
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Fig. 2. — The radial and vertical epicyclic frequencies for a neutron star rotating at a rate of 401 
Hz. The two thick dots indicate where the two epicyclic frequencies are separated by 1/2 the spin 
frequency (within grid errors). A torus whose circle of maximum pressure has radius rg = 11.27M 
(vertical line) would oscillate at frequencies 496 Hz and 694 Hz. The vertical dashed line indicates 
the radius (r = 8.72M) where the epicyclic frequencies are separated by the spin frequency. Also 
shown is ten times the Lense-Thirring frequency. The neutron star was modeled with equation of 
state FPS and has 1.22 solar mass. All frequencies are plotted as a function of the radius, given in 
units oiM{G/c^) = 1.22 x 1.5 km. 
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Fig. 3. — The frequencies (thick curves) of two g- modes of oscillation excited by parametric res- 
onance in a thin accretion disk. The two thick dots at r = 6.69M mark the boundary of the 
instability region — the resonant interaction with a rotating warp occurs to the right of the dots, 
where the lower frequency is less than the radial epicyclic frequency (dashed curve). The largest 
modulation of X-rays is expected close to the indicated boundary of the resonance region, i.e., at 
about 508 Hz and 737 Hz for the model computed (1.49 solar mass neutron star rotating at 406 
Hz, EOS AU). For a star of lower mass, i^iower and z^upper would have the values 500 and 700 Hz at 
a larger value of r/M. The computed frequencies of the two modes are uncertain by several percent 
because of the unknown rotation rate of the warp, which coincides at a poorly constrained radius 
with the Lense-Thirring frequency (also shown, in ten- fold magnification). 



